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I . INTRODUCTION 


1.1  Background 

The  original  EDIT-COMGEOM  code  was  written  by  Mr.  Larry  Bain  of 
Army  Materiel  Systems  Analysis  Agency:  the  version  of  the  EDIT-COMGEOM 

code  presented  in  this  report  is  an  improvement  of  the  original  code. 

1 2 

This  report  assumes  a familiarity  with  the  MAGIC  and  the  GIFT 
codes.  The  EDIT-COMGEOM  code  converts  target  description  data  used  as 
input  to  the  MAGIC  code  into  target  description  data  which  can  be  used 
as  input  to  the  GIFT  code. 

The  air  space  (01)  enclosing  the  target  used  by  the  MAGIC  code  can 
be  deleted  in  the  GIFT  code.  The  EDIT-COMGEOM  code  deletes  those  regions 
defined  as  "01"  air  spaces  and  when  possible  the  solids  associated  with 
those  regions  are  also  deleted.  The  GIFT  code  has  a new  solid  type 
called  RPP  which  requires  less  computer  storage  and  computer  run  time 
than  either  the  BOX  or  the  ARB  solid.  The  EDIT-COMGEOM  code  converts 
the  BOX  and  the  ARB  solids  which  are  equivalent  to  RPP's  into  RPP  solids. 

After  the  above  changes  are  made,  a renumbering  and  reordering  of 
the  target  description  data  for  the  MAGIC  code  is  required  before  the 
data  can  be  used  as  input  for  the  GIFT  code.  The  EDIT-COMGEOM  code 
does  the  renumbering  and  reordering. 

1.2  Contents  of  this  Report 

This  report  contains  three  sections  and  three  appendices.  The 
first  section  discusses  the  requirements  and  limitations  of  the  EDIT- 
COMGEOM  code.  The  second  section  describes  the  input  requirements  for 
the  code.  The  third  section  describes  the  output  of  the  code.  Appendix 
A is  a listing  (print-out)  of  the  code.  Appendix  B is  a listing  of  a 
sample  input.  Appendix  C is  a listing  of  the  output  from  the  sample 
input . 


NWC  Technical  Note  4565-3-71;  Volumes  I § II;  "MAGIC  Computer  Simula- 
tion, Volume  I,  User  Manual;  Volume  II,  Analyst  Manual,"  by  Armament 
Systems,  Inc.,  and  Propulsion  Development  Department,  May  1971. 
(Available  from  Department  of  the  Navy,  Naval  Weapons  Center,  China 
Lake,  California  93555). 

2 

Lawrence  W.  Bain  $ Mathew  J.  Reisinger,  "The  GIFT  Code  User  Manual; 
Volume  I Introduction  and  Input  Requirements,"  BRL  Report  No.  1802, 
April  1975. 
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II.  GENERAL  DISCUSSION 


2. 1 Required  Tapes 

The  EDIT-COMGEOM  code  requires  three  temporary  binary  tapes  on 
FORTRAN  units  10,  11,  and  12.  Disc  or  drum  may  be  substituted  for 
these  tapes.  The  tape  on  FORTRAN  unit  10  contains  the  region  table 
created  by  EDIT-COMGEOM  Code.  The  tape  on  FORTRAN  unit  11  contains 
the  MAGIC  Solid  Table,  and  the  tape  on  FORTRAN  unit  12  contains  the 
original  MAGIC  Region  Table. 

An  array  called  "IDENT"  is  used  to  store  the  Region  Identification 
Table.  It  is  dimensioned  for  1000  regions.  An  array  called  "ISOLID" 
flags  those  solids  to  be  kept  in  the  target  description  data.  It  is 
dimensioned  for  2500  solids.  The  Solid  Table,  Region  Table  and 
Region  Identification  Table  will  be  discussed  in  more  detail  in  later 
sections  of  this  report. 

2 . 2 Converting  the  Code  for  use  by  Different  Computers 

The  EDIT-COMGEOM  code,  listed  in  Appendix  A,  runs  on  either  the 
CDC  or  Univac  computers.  The  code  can  be  made  to  run  on  BRLESC 
(or  IBM)  computer  by  changing  the  comment  cards  labeled  in  card 
columns  73-80  with  "BRLESC"  C°r  "IBM")  to  executable  statements. 

Cards  labeled  in  columns  73-80  with  other  computers  must  be  removed  or 
changed  to  comment  cards. 

For  example,  the  arrowed  lines  in  the  listing  of  the  code  in 
Appendix  A are  as  follows: 


C 104  CALL  SET  CMC  (4)  BRLESC 

C 01 CODE  (80,32»S0L( 1)1  IREGt  ( I0PII1 v ISOLC 1 1 • 1*1 ,9 > BRLESC 

In  order  to  use  the  EDIT-COMGEOM  code  on  the  BRLESC  computer,  "C" 
would  be  removed.  The  four  lines  following  the  above  lines: 

104  WRITE  110,311  SOL  UNIVAC 

REWIND  13 

READ  110,40)  IREG,  CIOPI I ), ISOLCI1 , 1=1.9)  AND 

REWIND  10  CDC 


would  be  removed  or  changed  to  comment  cards. 
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III.  INPUT 


3.1  Control  Card 


The  control  card  is  the  first  data  card.  It  specifies  the  FORTRAN  input  and 
output  unit  numbers.  The  format  of  the  card  is  shown  in  Figure  1. 

3.2  MAGIC  Target  Description  Data 

The  MAGIC  target  description  data  follows  the  control  card.  The  form  of  the 
data  is  reported  in  reference  1,  "MAGIC  Computer  Simulation." 

The  data  consists  of: 

(a)  A Title  Card. 

(b)  A Target  Input  Constraints  Card. 

(c)  A set  of  solid  descriptions  called  a Solid  Table. 

(d)  A set  of  regions  described  in  terms  of  combinations  of  solids  called  a 
Region  Table. 

(e)  A set  of  alphanumeric  descriptions  of  a region  called  a Region  Identifica- 
tion Table. 

3.3  Target  Input  Constraints  Card 

The  input  data  for  the  EDIT-COMGEOM  code  uses  the  same  input  format  as  stated 
in  technical  note  "MAGIC  Computer  Simulation,"  except  for  the  Target  Input  Con- 
straints card.  The  format  of  the  Target  Input  Constraints  card  is  shown  in 
Figure  2. 

3.4  Order  of  Input 

A summary  of  the  order  of  the  data  for  EDIT-COMGEOM  code  is  shown  in  Figure  3. 

A listing  of  a sample  input  for  the  EDIT-COMGEOM  code  is  located  in  Appendix  B. 
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6-10 

11-15 

(BLANK) 

IN 

IOUT 

FORMAT:  (5X,  215) 

Name  Columns 

IN  6-10 

TOUT  11-15 


Figure  1.  Control  Card 


Remarks 

The  FORTRAN  unit  number  for  the  input. 

(If  "IN"  equals  0,  it  is  set  to  5 by  the 
code . ) 

The  FORTRAN  unit  number  for  the  output. 
(If  "IOUT”  equals  zero,  it  is  set  to  6 by 
the  code.) 


FORMAT:  (3110) 

Name  Columns 


NRPP 


NSOL 


NREGON 


1-10 


11-20 


21-30 


Remarks 

The  number  of  rectangular  parallelepipeds 
(RPPs)  used  to  describe  the  target's 
environment. 

The  number  of  solids  other  than  the  above 
RPPs  used  to  describe  the  target. 

The  number  of  regions  used  to  describe  the 
target . 


Figure  2.  The  Target  Input  Constraints  Card 
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£ 


£ 


£ 


L 


Z 


£ 


£ 


£ 


1/ 


Z 17 


Blank  Card  to 
Signal  End  of  Input 


Region  Identification  Table 


■ Blank  Card  (Optional) 
■1"  in  Columns  4 and  5 


Indicates  End  of  Region  Table 
Region  Table 


V 


Solid  Table 


Target  Constraints  Card 


Title  Card 


Control  Card 


Figure  3.  Card  Order  for  Input 


IV.  OUTPUT 


A listing  of  the  output  of  the  EDIT-COMGEOM  using  the  Sample  Input  of 
Appendix  B is  located  in  Appendix  C. 

4 . 1  Solid  Table 

A sample  line  of  output  of  the  Solid  Table  is: 

1RPR  -75,0000  75,0000  -36,0000  36,0000  12,0000  48,0000  2 3 

The  number  in  columns  1-3  (the  number  is  1 in  the  example  above)  is  the  new  solid 
number.  If  the  solid  number  is  greater  than  1000,  only  the  last  three  digits  are 
used.  The  number  in  columns  70-75  ("2"  in  the  example  above)  is  the  solid  number 
used  in  the  MAGIC  target  description  data.  The  number  in  columns  76-80  ("3"  in 
the  example  above)  is  the  card  or  line  count  of  the  GIFT  target  description  data. 


Those  ARB  or  BOX  solids  which  have  been  converted  to  an  RPP  have  the  following 
message  printed: 

BOX  3 HAS  BEEN  CONVERTED  TO  A RPR, 

4 . 2  Region  Table 


A sample  line  of  output  of  the  Region  Table  is: 

1 OR  1 -2  -6  -7  -ft  -9  -100R  30R  4 2 81 

The  number  in  columns  1-5  (the  number  is  1 in  the  example  above)  is  the  new  region 
number.  The  number  in  columns  71-75  ("2"  in  the  example  above)  is  the  region  number 
used  in  the  MAGIC  target  description  data.  The  number  in  columns  76-80  ("81”  in 
the  example  above)  is  the  card  or  line  count  of  the  GIFT  target  description  data. 

The  solid  numbers  in  columns  6-69  in  the  MAGIC  target  description  used  to 
describe  a region  are  replaced  with  the  new  solid  numbers  represented  in  the  new 
Solid  Table  of  the  output. 

4. 3 Region  Identification  Table 

A sample  line  of  output  for  the  Region  Identification  Table  is: 

1 too  0 BODY  2 96 

Columns  71-80  are  used  for  labeling  as  described  for  the  Region  Table  above, 

Any  alphanumeric  region  description  data  located  in  columns  71-80  of  the  MAGIC 
target  description  data  will  be  deleted. 

4.4  Summary  Table 

The  Summary  Table  contains  the  number  and  listing  of  solids  and  regions 
deleted  from  the  MAGIC  target  description  data.  In  the  sample  below, 1 solid  and  1 
region  were  deleted:  the  solid  deleted  and  region  deleted  were  both  numbered  "1": 
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NUMBER  OF  30LID8  DELETED  1 
1 

NUMBER  OF  RESIONS  DELETED  1 
1 
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Appendix  A:  A Listing  of  the  Code 


HAUTE R 

ED  IT-OOMGEOM  PROGRAM 

PROGRAM  DELETES  ALL  REGIONS  IDENTIFIED  AS  AIR  I 
AND  CHANGES  »NAGIC*  DOOIES  TO  'GIFT*  BOOIES 

ARRAYS  ARE  DIMENSIONED  FOR  1000  REGIONS  AND  2500  SOLIDS 

INPUT 

SPECIFY  INPUT  AND  OUTPUT  UNITS 

6-10  IN  - INPUT  UNIT  - IF  0 DEFAULTS  TO  CARD  READER  IN*5 
10-15  I OUT  - OUTPUT  UNIT  - IF  0 DEFAULTS  TO  PRINTER  10UT-6 
TITLE  CARD 
CONTROL  CARD 
SOL  10  TABLE 

REGION  TABLE  FOLLOWED  BY  -1  CARO 
I DENT  TABLE  FOLLOWED  BY  BLANK  CARD 

UNIT  IN  - INPUT  OLD  COM  GECM  DESCRIPTION 

UNIT  I OUT-  OUTPUT  NEW  COM  GEOM  DESCRIPTION  (EDITED) 

UNIT  10  - NEW  REGION  TABLE  (TEMPORARY  - BINARY) 

UNIT  11  - OLD  SOLID  TABLE  (TEMPORARY  - BINARY) 

UNIT  12  - OLD  REGION  TABLE  I TEMPORARY  - BINARY) 


DIMENSION  TITLE (18) 

DIMENSION  RPP»1B),  S0L(18),  1SOLID(2500) 

DIMENSION  I0P(9)t  IS0L(9) 

DIMENSION  KARD(ll),  IOENT 1 12. 1C00 ) 

EQUIVALENCE  !SOL( 2),S0LTYP) 

LOGICAL  DELETE*  PRNT 

DATA  BOX. ARB  /3HB0X* 3HARB  / BOX/ARB 

DATA  ATE  /4H8  / ARB 

DATA  ARS  /3HARS/ 

DATA  BLANK  / 3H  / 

DATA  END  / 3HEN0  / 

1 FORMAT! 5X  * 2 1 5 ) 

2 FORMAT  1 1H1.10HUN1T  IN  =>  .14/ 

1 1H  , 10HUNIT  OUT  *. 14// I 

10  FORMAT! 17A4.A2. 110) 

15  FORMATdH  ,17A4,AZ,U0) 

20  F0RMAT(3I10t40X, 110) 

25  FORMATdH  . 3110, 40X,  1101 

30  FORMAT (17A4.A2.2I5) 

31  FOBMATI I3.A3.16A4.2I5I 

32  FORMAT! 18, 9(A2, 15) ) BRLESC 

35  FORMATdH  ,17A4, A2.2I5) 

36  FORMATdH  , 1 3, A3,  16A4,  21 5 ) 

40  F0RMAT(I5,1X,9(A2,I5),1X,215) 

41  F0RMAT(5X,1X,9(A2,ISI,1X,2I5) 

42  FORMAT! 15, 70X, 15) 

43  F0RMAT(75X, 15) 

45  FORMATdH  , I 5»1X,9( A2, IS) ,1X, 2 15) 

46  FORMATdH  , 5X,  1X.9I  A2,  15 ) , IX,  21 5) 

47  FORMATdH  ,I5,70XVI5) 

48  FORMATdH  .75X.I5) 

50  F0RMAT(3I10.10X,7A4tA2,2I5) 

55  FORMATdH  ,3  I10,10X,7A4,  A2,  21 5) 
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INITIALIZE 

RB  AD  1 , IN  , I OUT 
IF  (IN.EQ.O)  IN  - 5 
IF  (IOUT.EQ.O)  (OUT  - 6 
WRITE (6. 2)  IN.IOUT 
NC ARO  - 0 
INCCRD  * l 

CALL  DISC  (10*2001  $ CALL  DISC  (11.2001  % CALL  OISC  (12.200) 
REWIND  10 
REWIND  11 
REWIND  12 

IF  (IN.NE.5)  REMIND  IN 
IF  (10UT.NE.6)  REWIND  (OUT 
PRNT-. FALSE. 

IF ( (OUT . EQ.6 ) PRNT*.TRUE. 

PRNT-.TRUE. 

CLEAR  IDENT  STORAGE 
DO  9D  1*1.12 
CO  90  L*  1 » 1000 

90  IDENT ( I.L )*0 

CLEAR  SOLID  TABLE 
CO  91  1*1,2500 

91  ISOL I D( I I *0 

READ  AND  WRITE  TITLE  CARD 

READ  (IN, 10 1 TITLE 

NCARO  * NC ARC  ♦ INCCRD 

IF(PRNT)  WRITE! IOUT, 15)  TITLE,  NCARD 

IF (.NOT.PRNT ) WRITE! IOUT, 10)  TITLE,  NCARD 

READ  CONTROL  CARD 

RB AO  ( I N,20 ) NRPP,  NSOL,  NREGON 
NS0L1D  = NRPP  ♦ NSOL 

READ  SOL  ICS  AND  STORE  ON  UNIT  11 

16  4 NRPP. EQ.O ) GOTO  101 
DO  100  1*1, NRPP 
READ  (IN, 30)  RPP 
WRITE  111)  RPP 

100  CONTINUE 

101  *S0L*0 

102  READ  (IN, 311  SOL 

IE  (SOLTVP. EC. BLANK)  GOTO  103 
NSOL-HSOL+l 

16  (NSOL. GT. NSOL)  GOTO  10* 

103  WRITE  (11)  SOL 
GOTO  102 

104  CALL  SET CWD  (4) 

DECODE  (80, 32, SOL(l))  IREG,  ( IOP( 1 ) , I SOL ( I ) * I *1 ,9) 

104  WRITE  (10,31)  SOL 
REWIND  10 

RBAD  (10,40)  IREG,  (I OP (1), ISOL (I), 1*1, 9) 

REWIND  10 
SOLTVP  - END 
WRITE  (11)  SOL 
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BRLESC 


HNl FSC 


bRLfcSC 

BRLESC 

UNI  VAC 
I BM 
AND 

coc 
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C READ  REGIONS  AND  STORE  ON  UNIT  12 

C 

110  WRITE  (12)  IRRG.I0P.IS0L 
IF  (IREG.LT.O)  GOTO  111 

READ  (IN, 40)  IREG.  ( I0P(I),  ISOLI I >, 1-1,9) 

GOTO  110 

C CHECK  FOR  REGION  RPP 

111  R8A0  (IN, 50)  KARO 
16  (KARDI 1I.NB.0)  GOTO  121 

REAO  IOENT  TABLE  ANO  STORE  IN  ARRAY  I0ENT 

120  READ  (IN, 50)  KARO 
IF  (KARO(l).EQ.OI  GOTO  200 

121  L*KARD(1) 

IDENTIltl 1*0 

IF  (KAROO).  EQ.ll  GOTO  120 
IF  ( K AROI 2 ) • EQ.999 ) GOTO  120 
00  1 22  1*1.11 
IOENT ( I ,L I - KAROII) 

122  CONTINUE 
IOENT (12,1 )x  KARO(l) 

GOTO  120 

REAO  OLD  REGION  TABLE  FROM  UNIT  12 
DELETE  AIR  1 

FLAG  SOLIOS  USED  IN  R6G10N  DESCRIPTION 
STORE  NEW  REGION  TABLE  ON  UNIT  10 

200  NBWR£G=0 
NREG*0 
REWINO  12 

210  REAO  (12)  IREG.IOP, ISOL 
IF  (IREG)  290,220*230 
SANE  REGION 

220  IF  (DELETE)  GOTO  210 
GOTO  250 

NEW  REGION 
230  NREG*NREG+1 

IF  ( IOENT (l.NREG).NE.O)  GOTC  240 
DELETE  * .TRUE. 

GOTO  210 

KEEP  THIS  REGION 
240  OELETE  * .FALSE. 

NEWREG  = NEWREG  * 1 
IOENT (1, NREG)  - NEWREG 
FLAG  SOLIOS  USEO 

250  00  251  1*1,9 
IF  (ISOLD). EQ.O)  GOTO  251 
L * I ABS  ( ISOL ( I ) ) 

ISOLIDIL)  * 1 

251  CONTINUE 
STORE  NEW  REGION 

WRITE  (10)  NEWREG, IOP, ISOL, IREG 
GOTO  210 

END  OF  REGIONS  FLAG 
290  WRITE  (10)  IREG, IOP, ISOL, IREG 

SET  UP  SOLIDS 
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NEWSOL  * 0 
00  300  1*1. NSOL  ID 
IF  (ISOLIO(I).EQ.O)  GOTO  300 
NEWSOL  * NEWSOL  ♦ 1 
I60LI0U1  * NEWSOL 

300  CONTINUE 
NBWRPP  - 0 

IF  <NRPP.EQ.0I  GOTO  310 
00  301  I*1,NRPP 

IR  (ISOLID! I). NE.O!  NEWRPP»NEWRPP+1 

301  CONTINUE 

WRITE  NEW  CONTROL  CARO 

310  NEWSOL  = NEWSOL  - NEWRPP 
NCARO  * NCARO  ♦ INCCRO 

IF(PRNT)  WRITE! I0UT»25)  NEWRPP, NEWSOL .NEWREG, NCARO 
IF ( .NOT.PRNT  I WRITE! I0UT.20!  NEWRPP , NEWSOL ,NEWREG .NCARO 

WRITE  NEW  RPP  CAROS 

REWINO  11 

IF  (NRPP.EQ.O)  GOTO  330 
CO  320  I RPP*1, NRPP 
READ  <111  RPP 

IF  <ISOLIO<  IRPPI.EQ.O)  GOTO  320 
NCARO  = NC ARC  ♦ INCCRO 
IF  <PRNT)  WRITE! I OUT, 35)  RPP, IRPP, NCARO 
IF (.NOT.PRNT!  WRITE! IOUT, 301  RPP, IRPP, NCARO 
320  CONTINUE 

WRITE  NEW  SOLID  TABLE 

330  N90L=NRPP 
340  READ  (11!  SOL 

IF  (SOLTYP.EG.ENDI  GOTO  400 
IF  (SOLTYP.NE. BLANK!  GOTO  350 
C SANE  SOLID 

IF  (DELETE!  GOTO  340 
GOTO  370 

C NEW  SOLI 0 

350  NSOL  * NSOL  ♦ 1 

IF  ( ISOLIO(NSOLI.NE.O!  GOTO  360 
DELETE  * .TRUE. 

GOTO  340 

C KEEP  THIS  SOLID 

360  DELETE  = .FALSE. 

NEWSOL  * M00(  ISOLID(NSOL)  , 1000  ! 

C CHANGE  'MAGIC*  ARS  TO  *GIFT*  ARS 

IF (SOLTYP.NE. ARSI  GOTO  364 
DO  362  1*4,18 

If (SOL! I I.NE. BLANK  I GOTO  370 
362  CONTINUE 

READ! 11)  SOL 
SOLTYP-ARS 
GOTO  370 

C CHECK  BOX  OR  ARO  FOR  OBGENERATION  TO  AN  RPP  BODY 

364  IF (SOLTYP.NE. BOX)  GOTO  369 

CALL  CKBOX  < SOL .NEWSOL , NCARO* INCCRO, NSOL, IOUT, PRNT! 

GOTO  340 


16 


non  o noo  o o o o o 


3 


36$  IF (SOLTYP.NE. ARB)  60T0  370 

in  I SOL (3 ) .NE. ATE.AN0.S0LI3) .NE.BLAMK)  GOTO  370 
CALL  CKARB  I SOL, NEWSOL, NCARO, INCCRO, NSOL, I OUT, PRNT) 

GO  70  340 

C WRITE  NEW  SOL  10 

370  NCARO  * NCARO  ♦ INCCRO 

IFIPRNT)  WRITEI I0UT,36I  NEWSOL, ISOL 1 1 ), 1*2,18) , NSOL, NCARO 
Ifll. NOT.  PRNT)  WRITEI  I0UT,31)  NEWSOL.  I SOLID  , 1-2  * IB)  ,NSOL, NCARO 
GOTO  340 

WRITE  NEW  REGION  TABLE 
400  RGW1N0  13 

410  RE AO  (101  NEWREG, IOP, I SOL, IREG 
IF  (NEWREG.LT. 0)  GOTO  490 
00  420  1=1.9 

IF  ( ISOL( Il.EQ.O)  GOTO  420 
L * I ABS  (ISOL (I)) 

ISOLII)  = ISIGN  ( ISOLIOIL),  (SOLUM 
420  CONTINUE 

NCARO  * NCARO  ♦ INCCRO 
IF  IIREG.EQ.O)  GOTO  43G 
NREG  - IREG 

IF ( PRNT I WRITEI IOUT, 45)  NEWREG, ( IOP (I) , I SOL ( I) . I = 1 ,9 > , NREG, NCARO 
IF  I. NOT. PRNT)  WRITEI I OUT, 40 INEWREG, ( I OP ( 1 1 , 1 SOL ( I ) , 1 = 1 ,9) .NREG, 

1 NCARO 

GOTO  410 

430  IFIPRNT)  WRITE! IOUT, 46)( 10P( I ), ISOLII ), 1=1, 9), NREG, NCARO 

IF (.NOT. PRNT I WRITEI IOUT, 46 M I OP II), ISOL (11,1  = 1,9), NREG, NCARO 
GOTO  410 

END  OF  REGION  FLAG 
490  NCARO  = NCARO  ♦ INCCRO 

IFIPRNT)  WRITEI I0UT.47)  NEWREG, NCARO 
IF  I . NOT. PRNT)  WRITEI IOUT, 42)  NEWREG, NCARO 
BLANK  CARO  FOR  REGION  RPP 
NCARO  = NCARO  ♦ INCCRO 
IFIPRNT)  WRITEI IOUT, 48)  NCARO 
IF ( . NOT. PRNT)  WRITEI IOUT, 43)  NCARO 

WRITE  NEW  IOENT  TABLE 

00  $00  L=1 , NREGON 
IF  I IOENT 1 1, L I .EQ.O ) GOTO  500 
NCARO  * NCARO  ♦ INCCRO 

IFIPRNT)  WRITEI IOUT, 55)  I IOENT! I, L), 1=1, 12), NCARO 
If (.NOT. PRNT)  WRITEI IOUT, 50  1 1 IOENT! I, L), 1 = 1, 12) .NCARO 
500  CONTINUE 

END  OF  TABLE  FLAG 
NCARO  * NCARO  ♦ INCCRO 
IFIPRNT)  WRITEI IOUT, 48)  NCARO 
IF{. NOT. PRNT)  WRITEI (OUT, 43)  NCARO 

LIST  SOLIOS  DELETED 

NSOL  - 0 

00  1000  I*1,NS0LI0 
If  IISOLIO(I).NE.O)  GOTO  1000 
NSOL  - NSOL  ♦ 1 
ISOLIOINSOL)  =i  I 
1000  CONTINUE 
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PRINT  1001  , NSOL 

IF  INSOL. NE. 01  PRINT  1002  . < IS0L 101 11.1=1 t NSOL I 

1001  FORMAT ( 1H1, ' NUMBER  OF  SOLIDS  DELETE0' ,15  I 

1002  FORMAT ( 1H  ,10151 


LIST  REGIONS  DELETED 


NR  EG  * 0 

DO  2000  1*1 , NREGON 

IF  1 1 DENT  11*11. NE .01  GOTO  2000 

NREG  * NR EG  ♦ 1 

I&OL I D( NREG 1 * I 

CONTINUE 

PRINT  2001  , NREG 

PRINT  2002  , I I SOL  1 01  I ), 1*1, NREG) 


2001  FORMAT ( 1HO, * NUMBER  OF  REGIONS  OELETED  * , 15  ) 

2002  FORMATdH  ,1015) 


STOP 

END 


SUBROUTINE  CKARBt SOL,NEWSOL,NCARD, INCCRO.NSOL , I OUT , PRNT ) 

THIS  SUBROUTINE  CHANGES  AN  ARB  TO  A RPP 
DOUBLE  PRECISION  X,Y,Z,A,B,C 

DIMENSION  S0L(1B),X(8),Y(8),Z(8),R(6),IV(4,6),A(6) ,B(6> ,C(6) 

DIMENSION  TEMPdCI  UNIV/COC 

LOGICAL  PRNT 
DATA  ARB  /3HARB/ 

DATA  ATE  /4H8  / 

DATA  BLANK  /4H  / 

DATA  RPP  /3HRPP/ 

DATA  IV  / 1,2, 3, A,  5, 6, 7, 8,  1,5,8, A,  2, 3,7,6,  1,2, 6, 5.  4, 3, 7, 8/ 

1 F0RMATI2A3.16A4) 

2 FORMAT ( 6F10 • 0 1 

3 F0RMATd0X,6F10.0) 

4 FORMATdH  , 13,  A3,  A4.6F10.4, 2151 

5 FORMATdH  , 1 3 , 7X, 6F 10. A,  2 15  1 

6 FORMAT ( 5H  ARB  , I3.29H  HAS  BEEN  CONVERTED  TO  A RPP.) 

7 FO RMATI6 (1X,4 1 1 ) ) 

8 FORMAT!  10X,  6 dX,  AI 1 1 1 

9 FORMAT ( 15A4) 

10  FORMATdH  , 1 3,7X, 6!  IX, 41 1 1 , 30X, 21 5 1 

11  FORMAT! I3,A3,A4,6F10>4,2I5) 

12  FORMAT! I3,7X,6F10.A,2I5) 

13  F0RMAT(13,7X,6!1X,AI1),30X,2I5) 

T0L**O001 

HS0L=SOL 13) 

CONVERT  TO  FLOATING  POINT  FORMAT 


REWIND  12  IBM 

DO  100  1*1, 8, 2 

EM CODE (60 ,9 , TEMPI  1 ) ) I SOL !J1»J*4»18)  UNI  VAC 

DBC0DE(60 ,2 > TEMP! 1 ) I X! 1 ) ,V! I ) ,Z( I ) ,X( 1*1 ) , V! 1*1 ) ,2  ! I d I ANO  COC 

DBCOOE!  60, 2,  SOL!  A I I X(I),Yd),Zd),X(I*lt,Y(I+l),ZII  + l>  BRLESC 

WR  1TE!  12 ,1 ) SOL  IBM 

100  IF(N*NE.A.OR.HSOL.NE.ATE)  READ! 11)  SOL 
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IMIHSOL.EQ.BLANK)  WRITE! 12,1)  SOL  IBM 

RBWINO  12  IBM 

DO  101  I-1,B,2  IBM 

101  REA0I12,3)  xmiYIIM(IM(Hl|,Y(M),Z(mi  IBM 

IFIHSOL.EQ.BLBNK)  READ! 12, B)  IV  IBM 

IFIHSOL.EQ.BLANK)  DECOOE! 30, 7, SOL !4 I)  IV  BRLESC 

IFIHSOL.EQ.ATEI  GOTO  110  UNIVAC 

ENCODE 1 3 2,9 tTEMP| 1))( SOLI J I , J*4, 11 1 

OECODEIBQfT, TEMPI  II ) IV  CDC 


C FINO  EQUATION  OF  FACES,  TEST  FOR  RPP 

110  NZ-0 

DO  120  1 = 1,6 

CALL  EQFACE  I AP.BP, CP, XI L I , YIL I ,Z!L> , 1 ,2, 31 
IFlABSflAim.LE.TOLI  NZ-NZ+1 
IFIABSiail  H.LE.TOL)  NZ»NZ+l 
120  IFIABSICIIII.LE.TOL)  NZ-NZ+1 
IRINZ.EQ.12)  GOTO  200 
C FAILED,  REMAINEO  AN  ARB  1 

NCARD-NCARO+INCCRO 

IF! PANT  I WRITE! I0UT,4 INEWSOL, ARB.HSOL, IX!  1 1 , Y!II , Z (II , 1*1 ,2 > ,NSOL, 

1 NCARD 

IF ! « NOT. PRNT  > WRITE! I OUT ,11 INEWSOL, ARB ,HSOL ,(XII),Y(I),Z(I),I*1,2) 
1 ,NSOL, NCARD 

DO  130  1=3, 8, 2 
NC ARD-NCARD+ INCCRD 

IF ( PRNT I WRITE! IOUT, 5 ) NEWSOL ,X ( I),Y(1),ZII),X(I+1),Y|I+1),Z(I+1), 

1 NSOL, NCARD 

IF ( .NOT . PRNT )WR ITE! IOUT ,12)NEWS0L,X!I ) , Y( I ) ,Z< I)  »X ! 1 + 1 ) ,Y!I+1>, 

1 Zl I+l), NSOL, NCARD 

130  CONTINUE 

IFIHSOL.EQ.ATEI  RETURN 
NCARD-NCARD+ INCCRD 

IF! PRNT I WRITE! IOUT, 101  NEWSOL, IV«NSOL, NCARD 

IF ( .NOT. PRNT I WRITE! IOUT, 131  NEWSOL, IV, NSOL, NCARD 

RETURN 

C CONVERT  TO  RPP 

200  DO  210  1*2,8 

IFIXi II.NE.XIll )X!2)-X< I 1 
IFIY!I).NE.YI1|)YI2)«Y|I| 

210  IF(ZII).NE.Z(1))Z(2)=ZII) 

RI1)«ANIN1IXI1),X|2I) 

RI2)=AHAX1IXI1),XI2I) 

RI3I-AMIN1IYI 1),YI 2)1 
R44)»AMAX1IYI1),Y!2) I 
R15)-AMIN1!ZI1),Z!2)I 
R1 6 )• AMAX1I Z!1I,ZI2I! 

NCARD-NCARO+ INCCRD 

IF ( PRNT I WRITE! IOUT, 41  NEWSOL, RPP, BLANK, R, NSOL, NCARD 

IF  I .NOT. PRNT I WRITEIIOUT.llI  NEWSOL.RPP, BLANK, R, NSOL, NCARD 

WRITEI6.6I  NSOL 

RETURN 

END 


SUBROUTINE  CKBOX  {SOL, NEWSOL, NCARO, INCCRD, NSOL, IOUT, PRNT) 
C 

C THIS  SUBROUTINE  CONVERTS  A BOX  TO  A RPP 
C 

DIMENSION  SOLI  18 1 ,BI 12 1, TEMPI 10I,RI 6) 

LOGICAL  PRNT 
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DATA  BOX  /3HB0X/ 

MTA  RPP  /3HRPP/ 

1 F0RHAT12 A3,16A4) 

2 FORMATI6F10.0) 

3 FORMAT! 10X,6F10.0J 

4 FORMAT! 1H  , I3»A3,4X,6F10.4,2I5) 

5 FORMAT ! 1H  , 13. 7X , 6F10.4, 2 15 1 

6 FORMAT! 5H  BOX  , I3,29H  HAS  BEEN  CONVERTED  TO  A RPP.) 

7 FORMAT! 15A4) 

B F0RNATII3,A3,4X,6Fl0. 4,215) 

9 FORMAT! 13, 7X.6F10. 4,215) 

C CONVERT  TO  FLOATING  POINT  FORMAT 


ENCODE!  60, 7,  TEMPI  L) ) I SOL!  I ), 1=4,18)  UNI  VAC 

DECODE  1 60, 2,  TEMP!  1))(B(I),I=1,6)  AND  CDC 

DECODE!  6 0,2,  SOU  4)  )!B(  I ),  1-1,6)  BRLESC 

RBHIND  12  IBM 

MR  I T E ( 12 • 1 ) SOL  IBM 

RBAD(Ll)  SOL 

ENCODE  160,7 ,TEMP ! 1 ) ) ( SOLI  I ), 1 = 4,18)  UN  I VAC 

DEC0DE(60,2, TEMP! 1))(B(I),1=7,12)  AND  CDC 

DECODE! 60, 2, SOL! 4)  MB!  I),  1=7,12)  BRLESC 

UR  I T E ( 12 , 1 ) SOL  IBM 

REMIND  12  IBM 

READ (12,3)  B IBM 


CHECK  FOR  RPP 

NZ-0 

DO  100  1=4,12 

100  IP ( B ( I ) . EQ.O .0 ) NZ  =NZ ♦ 1 
IFINZ.EQ.6)  GOTO  200 

FAILED,  REMAINS  A BOX 
NCARD=NCARD+  INCCRD 

IF(PRNT)  MRITE! I OUT, 4)  NEMSOL , BOX, ( 8 ( I ) , I =1 ,6 ) , NSCL ,NCA RD 

IF ( .NOT.PRNT I MRITE! I0UT,8)  NEMSOL, BOX, IB! I ) ,1=1,6) ,NSOL,NCARD 

NC  ARD*NC ARD+ INCCRD 

IF(PRNT)  MRITE! I0UT.5)  NEMSCL , I Bl I) • I =7, 12 ) ,NSOl ,NCARD 

IF (.NOT.PRNT)  MRITE!  I0UT.9)  NEMSOL , I B ( 1 1 , I =7, 12 ) ,NSOL ,NCARD 

RBTURN 

CONVERT  BOX  TO  RPP 
200  B(4)=B(1)*B!4)*BI7)*BI10) 

Bi5)=B(2)*B(5l+B(81*B(  111 
B(6)=B(3l+B(6)+B(9)+8( 121 
R ( 1 1 = AM  I N 1 ! B ( 1 ) , B ( 4 ) ) 

R12)«AMAX1(B( 11,8(41 ) 

R! 3 ) = AMI N1 I B ( 2 ) • B( 5 ) ) 

R(4)=AMAX1(B(2),B(5) I 
R15)*AMIN1(B(3),B(6) ) 

R(6)=AMAXl(BI3),B(6l I 
NCARD=NCARD+ INCCRD 

IFIPRNT)  MRITE! I0UT.4)  NEMSOL , RPP, R,NSOL,NCARD 

IF!. NOT.PRNT)  MRITE! IOUT, 8 I NEMSOL, RPP, R.NSOL ,NCARD 

MR ITE ( 6, 6 ) NSOL 

RETURN 

END 


SUBROUTINE  EQFACE  I X, Y, Z , A, B, C,K, L ,M I 

THIS  SUBROUTINE  CALCULATES  A,B,C  OF  THE  EQUATION  OF  FACE  AX*BY+CZ=1.0 
If  ALSO  CALCULATES  THE  POINT  OF  INTERSECT  OF  THREE  PLANES 
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DOUBLE  PRECISION  A,8,C»D,X,Y,Z 
DIMENSION  XI8>,Y(8),*I8J 

D-.X(K)*Y(L>*ZtH>+XIM>*Y(K»*Z<L>*XIU*YtMI^ZIK»-X«M)*Y(U*Zm-X(L> 
1 *Y  ( K ) *Z  ( N I — X ( K I *Y  ( M I *Z  I L I 
IF  ( D.  EQ.'J.O  ) D». OOOOl 

A>  mK)*(Z(L)-Z(MII«YfU*(Z(M>-ZtK))-frY(H)*(Z<K)-Z(Lm/D 

B^(Z(KIMX(LI-X(M|)«mi*(X(X|>X(K)|tZ(M|*|X(K|'XlLl)l/D 

C^(X(KI*IY(L|-YIM| )*X(LI*(Y(M |*Y IK ) ) ♦X ( M ) * ( Y ( K l-Y(L) I I/O 

RETURN 

END 
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Appendix  B:  A Sample  Input  for  Code 


I/O  5 6 

SAMPLE  INPUT  USED  IN  NUC  DOCUMENTATION 


l 24  12 


-10000. 

10000. 

-10000. 

10000. 

-10000. 

10000. 

2 BOX 

75. 

-36. 

12. 

-150. 

0. 

0. 

BODY 

2 

0. 

72. 

0. 

0. 

0. 

36. 

3 BOX 

74. 

-35. 

13. 

-148. 

0. 

0. 

(1.0) 

3 

0. 

70. 

0. 

0. 

0. 

34. 

4 ARB 

75. 

-36. 

12. 

75. 

36. 

12. 

FRONT 

4 

75. 

36. 

48. 

75. 

-36. 

48. 

4 

100. 

0. 

12. 

100. 

0. 

12. 

4 

100. 

0. 

12. 

100. 

0. 

12. 

4 

1234  6435 

6128  6237 

7415 

7415 

5 ARB 

-75. 

-36. 

12. 

-75. 

36. 

12. 

REAR 

5 

-75. 

36. 

48. 

-75. 

-36. 

48. 

5 

-100. 

-24. 

12. 

-100. 

24. 

12. 

5 

-100. 

24. 

20. 

-100. 

-24. 

20. 

5 

1234 

5678  3487 

1265  2376 

1485 

6 ELL 

20. 

0. 

48. 

• 

0 

CM 

1 

0. 

48. 

8UBBLE 

6 

50. 

7 ELL 

7 0. 

0. 

48. 

24. 

0. 

C. 

(1.0) 

7 

14. 

8 RCC 

60. 

-36. 

12. 

-0. 

8# 

0. 

wheel 

8 

12. 

9 RCC 

60. 

36. 

12. 

0. 

-8. 

0. 

WHEEL 

9 

12. 

10  RCC 

—60. 

-36. 

12. 

0. 

8. 

0. 

WHEEL 

10 

12. 

11  RCC 

-60. 

36. 

12. 

c. 

-8. 

u ■ 

WHEEL 

11 

12. 

12  BOX 

-70. 

-20. 

15. 

40. 

0. 

c. 

ENGINE 

12 

0. 

40. 

0. 

0. 

0. 

3 0. 

13  RAW 

-70. 

-20. 

45. 

0. 

0. 

-10. 

(ENGINE) 

13 

0. 

10. 

0. 

4C. 

0. 

0. 

14  RAW 

-70. 

20. 

45. 

0. 

0. 

-10. 

(ENGINE) 

14 

0. 

-10. 

0. 

40. 

0. 

c. 

15  ARB 

-70. 

-10. 

45. 

-70. 

10. 

45. 

(ENGINE) 

15 

-70. 

0. 

35. 

-7C. 

0. 

35. 

15 

-30. 

-10. 

45. 

-30. 

10. 

45. 

15 

-30. 

0. 

35. 

-30. 

0. 

35. 

15 

3124 

7658  1375 

2376  1265 

1265 

16  ARS 

16 

4 

5 

16 

-70. 

-20. 

15. 

-70. 

-20. 

15. 

1 

16 

-70. 

-20. 

15. 

-70. 

-20. 

15. 

2 

16 

-70. 

-20. 

15. 

3 

16 

-70. 

-20. 

15. 

-70. 

-10. 

15. 

4 

16 

-70. 

-10. 

25. 

-70. 

-20. 

35. 

5 

16 

-70. 

-20. 

15. 

6 

16 

-30. 

-20. 

15. 

-30. 

-10. 

15. 

7 

16 

-30. 

-10. 

25. 

-30. 

-20. 

35. 

8 

16 

-30. 

-20. 

15. 

9 

16 

-30. 

-20. 

15. 

-30. 

-20. 

15. 

10 

16 

-30. 

-20. 

15. 

-30. 

-20. 

15. 

11 

16 

-30. 

-20. 

15. 

12 

17  ARS 

17 

5 

4 

17 

-70. 

20. 

IS. 

-70. 

20. 

15. 

1 

17 

-30. 

20. 

15. 

-30. 

20. 

15. 

2 

23 


2 


17 

-70. 

20. 

15. 

-70. 

10. 

15. 

17 

-30. 

10. 

15. 

-30. 

20. 

15. 

17 

-70. 

20. 

15. 

-70. 

10. 

25. 

17 

-30. 

10. 

25. 

-30. 

20. 

15. 

17 

-70. 

20. 

15. 

-70. 

20. 

35. 

17 

-30. 

20. 

35. 

-30. 

20. 

15. 

17 

-70. 

20. 

15. 

-70. 

20. 

15. 

17 

-30. 

20. 

15. 

-30. 

20. 

15. 

18 

REC 

0. 

0. 

24. 

0. 

C. 

28. 

TRUNK 

18 

0. 

7.5 

0. 

5. 

0. 

'J  m 

19 

SPH 

0. 

0. 

52. 

5. 

HEAD 

20 

TEC 

0. 

-7.5 

49. 

20. 

0. 

-12. 

arm 

20 

0. 

0. 

3. 

0. 

2. 

0. 

20 

2. 

21 

TEC 

0. 

7.5 

49. 

20. 

0. 

-12. 

ARM 

21 

0. 

0. 

3. 

0. 

2. 

0. 

21 

2. 

22 

TRC 

-2. 

-4.5 

27. 

32. 

0. 

-12. 

LEG 

22 

3. 

2. 

23 

TRC 

-2. 

4.5 

27. 

32. 

0. 

-12. 

LEG 

23 

3. 

2. 

24 

TOR 

21.5 

C. 

37. 

1. 

c. 

Cm 

STEERI NG 

24 

8. 

1. 

khe  el 

25 

ARB 

21.5 

-6. 

33.5 

21.5 

6. 

33.5 

CENTER 

25 

21.5 

0. 

44. 

41. 

3. 

37. 

STEERING 

25 

21.5 

-6. 

33.5 

21.5 

6. 

33.5 

WHEEL 

25 

21.5 

0. 

44. 

4C. 

0. 

37. 

25 

3127 

2146  4328 

1345 

3127 

3127 

1 

1 

-2 

-4 

-5 

-6  - 0 

-9 

-Ij 

-11 

2 OR 

2 

-3 

-7 

-9 

-9  -1C 

- 1 ICR 

4^R 

5 

3 

6 

-7 

-2 

4 

8 

5 

9 

6 

10 

7 

11 

8 OR 

3 

-18 

19 

-2G 

-21  -22 

-23 

-2  4 

-25 

-8 

-9 

10 

-li 

-12CR  7 

-18 

-19 

-20 

-21 

-24 

250R 

130R 

14GR  150R 

16CR 

17 

9 

3 

12 

13 

-14 

-15  -16 

-17 

10  OR 

180R 

190R 

200R 

210R 

220R  23 

11 

3 

24 

25 

12 

3 

25 

-1 

1 

Cl 

OUTSIDE  i 

AIR 

2 

100 

BODY 

3 

101 

BUBBLE 

4 

151 

WHEEL  RIGHT 

FRONT 

5 

152 

WHEEL  LEFT 

FRONT 

6 

153 

WHEEL  RIGHI 

REAR 

7 

154 

WHEEL  LEFT 

REAR 

8 

02 

AIR  INSIOE 

9 

200 
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